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Abstract 

By use of a stochastic generalization of the Hodgkin-Huxley model we 
investigate both the phenomena of stochastic resonance (SR) and coher- 
ence resonance (CR) in variable size patches of an excitable cell mem- 
brane. Our focus is on the challenge how internal noise stemming from 
individual ion channels does affect collective properties of the whole en- 
semble. We investigate both an unperturbed situation with no applied 
stimuli and one in which the membrane is stimulated externally by a pe- 
riodic signal and additional external noise. For the nondriven case, we 
demonstrate the existence of an optimal size of the membrane patch for 
which the internal noise causes a most regular spike activity. This phe- 
nomenon shall be termed intrinsic CR. In presence of an applied periodic 
stimulus we demonstrate that the signal-to-noise ratio (SNR) exhibits SR 
vs. decreasing patch size, or vs. increasing internal noise strength, respec- 
tively. Moreover, we demonstrate that conventional SR vs. the external 
noise intensity occurs only for sufficiently large membrane patches, when 
the intensity of internal noise is below its optimal level. Thus, biological 
SR seemingly is rooted in the collective properties of large ion channel 
ensembles rather than in the individual stochastic dynamics of single ion 
channels. 

PACS number(s): 87.10.+e, 87.16.-b, 05.40.-a 

During the last decade, the effect of stochastic resonance (SR) - a coopera- 
tive phenomenon wherein the addition of external noise improves the detection 
and transduction of signals in nonlinear systems (for a comprehensive survey 
and prominent references, see Ref. [ID) - has been studied experimentally and 
theoretically in various biological systems ||, [| @! §|> P|- For example, SR has 
been experimentally demonstrated within the mechanoreceptive system in cray- 
fish ||, in the cricket cereal sensory system [|J, for human tactile sensation Q, 
visual perception ||, and response behavior of the arterial baroreflex system 
of humans Q . The importance of this SR-phenomenon for sensory biology is 
by now well established; yet, it is presently not known to which minimal level 
of the biological organization the stochastic resonance effect can ultimately be 
traced down. Presumably, SR has its origin in the stochastic properties of the 



ion channel clusters located in a receptor cell membrane. Indeed, for an artificial 
model system Bezrukov and Vodyanoy have demonstrated experimentally that 
a large parallel ensemble of the alamethicin ion channels does exhibit stochastic 
resonance Q. This in turn provokes the question whether a single ion channel 
is able to exhibit SR, or whether stochastic resonance is the result of a collective 
response from a finite assembly of channels. Stochastic resonance in single, bio- 
logical potassium ion channels has also been investigated both theoretically || 
and experimentally |9| . This very experiment did not convincingly exhibit SR in 
single voltage-sensitive ion channels. Nevertheless, the SR phenomenon can oc- 
cur in a single ion channel if only the parameters are within a regime where the 
channel is predominantly dwelled in the closed state, as demonstrated within a 
theoretical modeling for a potassium Shaker channel This prominent result, 
i. e. the manifestation of SR on the single- molecular level, is not only of academic 
interest, but is relevant also for potential nanotechnological applications, such as 
the design of single- molecular biosensors. The origin and biological relevance of 
SR in single ion channels, however, remains still open. Indeed, biological SR can 
be a manifestation of collective properties of large assemblies of ion channels of 
different sorts. To display the phenomenon of excitability these assemblies must 
contain a collection of ion channels of at least two different kinds - such as, e.g., 
potassium- and sodium-channels. The corresponding mean-field type model has 
been put forward by Hodgkin and Huxley in 1952 |l(J by neglecting the meso- 
scopic fluctuations which originate from the stochastic opening and closing of 
channels. SR due to external noise in this primary model and related models of 
excitable dynamics has extensively been addressed jll|]. These models further 
display another interesting effect in the presence of noise, namely so termed co- 
herence resonance (CR) even in absence of an external periodic signal the 
stochastic dynamics exhibits a surprisingly more regular behavior solely due to 
an optimally applied external noise intensity. A challenge though still remains: 
Does internal noise play a constructive role for SR and CR? Internal noise is 
produced by fluctuations of individual channels within the assembly, and dimin- 
ishes with increasing number. For a large, macroscopic number of channels this 
noise becomes negligible. Under the realistic biological conditions, however, it 
may play a crucial role. 

Our starting point is due to the model of Hodgkin and Huxley jl3|, i.e. 
the ion current across the biological membrane is carried mainly by the motion 
of sodium, Na + , and potassium, K+, ions through selective and voltage-gated 
ion channels embedded across the membrane. Besides, there is also a leakage 
current present which is induced by chloride and remaining other ions. The ion 
channels are formed by special membrane proteins which undergo spontaneous, 
but voltage-sensitive conformational transitions between open and closed states 
[ pi] . Moreover, the conductance of the membrane is directly proportional to 
the number of the open ion channels. This number depends on the potential 
difference across the membrane, V. The different concentrations of the ions 
inside and outside the cell are encoded by corresponding reversal potentials 
E^a — 50mV, Ek — — 77mV and E^ — — 54.4mV, respectively, which give the 
voltage values at which the direction of partial ion currents is reversed fll3|. 
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Taking into account that the membrane possesses a capacitance C = 1/zF/cm 2 , 
Kirchhoff's first law reads in presence of an external current I ex t(t) stimulus: 

C i v + Gk W (Y - £ k) + G Na (m, h) {V - £ Na ) + G L (V- E L ) = I cxt {t) . 
at 

(1) 

Here, G^ a .{m, h), Gk(ti) and Gl denote the conductances of sodium, potassium, 
and the remaining other ion channels, respectively. The leakage conductance 
is assumed to be constant, Gl = 0.3mS/cm 2 ; in contrast, those of sodium and 
potassium depend on the probability to find the ion channels in their open con- 
formation. To explain the experimental data, Hodgkin and Huxley did assume 
that the conductance of a potassium channel is gated by four independent and 
identical gates. Thus, if n is the probability of one gate to be open, the probabil- 
ity of the K + channel to stay open is Pk — n • Moreover, the gating dynamics 
of sodium channel is assumed to be governed by three independent, identical 
gates with opening probability m and an additional one, being different, pos- 
sessing the opening probability h. Accordingly, the opening probability of Na + 
channel (or the fraction of open channels) reads PNa = rv?h. The conductances 
for potassium and sodium thus read 



G K (n)= 5 r x ™ 4 , G^{m 1 h) = g^ m i h 1 (2) 

where g^& x = 120mS/cm 2 and g^ ax = 36mS/cm 2 are the maximal conduc- 
tances. Furthermore, the gating variables (probabilities) m, h and n obey the 
two-state, "opening-closing" dynamics, 

m = a m (V) (1 - m) - m (V) m , 
h = a h (V) (l-h)-/3 h (V) h, (3) 
h = a n (V) (l-n)-(3 n (V) n, 

with the experimentally determined voltage-dependent transition rates, reading 

""•< y > = l-expRV^'o)/^ ' " 4 »"|-<V + 65)/18] , 

a h (V) = 0.07 exp[-(V + 65)/20], (3 h {V) = {l + exp[-(V + 35V10]}- 1 , 

(4) 

^ = MV) = °- 125 cxph(F + 65)/80] • 

These nonlinear Hodgkin-Huxley equations ([[])- (f|) present a cornerstone model 
in neurophysiology. Within the same line of reasoning this model can be gen- 
eralized to a mixture of different ion channels with various gating properties 
[Et| E5l. An essential drawback of this model, however, is that it operates with 
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the average number of open channels, thereby disregarding corresponding num- 
ber fluctuations (or, the so-called channel noise |l6|]). Thus, it can strictly be 
valid only within the limit of very large system size. We emphasize, however, 
that the size of an excitable membrane patch within a receptor cell is realis- 
tically finite. As a consequence, the role of internal fluctuations cannot be a 
priori neglected; as a matter of fact, as shown below, they do play a key role 
for SR and CR. 

The role of channel noise for the neuron firing has been studied by Lecar and 
Nossal as early as in 1971 @. The corresponding stochastic generalizations of 
Hodgkin- Huxley model (within a kinetic model which corresponds to the above 
given description) has been put forward by DeFelice et al. jl8| and others; see 
[ fl6| for a review and further references therein. The main conclusion of these 
previous studies is that the channel noise can be functionally important for 
neuron dynamics. It particular, it has been demonstrated that channel noise 
alone can give rise to a spiking activity even in the absence of any stimulus 

The stochastic kinetic scheme |16|, uM, however, necessitates extensive nu- 
merical simulations . To aim at a less cumbersome numerical scheme we use 
a short-cut procedure that starts from Eq. (|J) in order to derive a corresponding 
set of Langevin equations for a stochastic generalization of the Hodgkin-Huxley 
equations of the type put forward by Fox and Lu Q . Following their reasoning 
we substitute the equations (0) with the corresponding Langevin generalization: 



m = a m (V) (1 - m) - (3 m (V) m + £ TO (t) , 

h = a h (V){l-h)-0h(V)h + t h (t) t (5) 
n = a n (V) (1 - n) - (3 n (V) n + £„(*) , 

with independent Gaussian white noise sources of vanishing mean. The noise 
autocorrelation functions depend on the stochastic voltage and the correspond- 
ing total number of ion channels as follows: 



A*Na \a m + Pm ) 

(UW)) = ^- . ah f" . S(t - t>) , (6) 
A^Na (ah + Ph ) 

N K (a n + Pn) 

In order to confine the conductances between the physically allowed values be- 
tween (all channels closed) and g max (all channels open) we have implemented 
numerically the constraint of reflecting boundaries so that m(t), h(t) and n(t) 
are always located between zero and one [p(| . 

Moreover, the numbers iVN a and Nk depend on the actual area S of the 
membrane patch. With the assumption of homogeneous ion channels densities, 
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PNa = 60^m 2 and = 18y«m 2 , the following scaling behavior follows: 

A'Na — PN&S, N K = PK S. (7) 

Upon decreasing the system size S, the fluctuations and, hence, the internal 
noise increases. 

Before integrating the system of stochastic equations (0), (j|), (^) numeri- 
cally, the external stimulus I e xt(t) must be specified. We take a periodic stimulus 
of the form 



I ext (t) =A sin(m) + r?(t) , (8) 

where the sinusoidal signal with amplitude A and frequency Q is contaminated 
by the Gaussian white noise rj{t) with the autocorrelation function 

(T?(*)TJ(0> = 2 ^ext *(i - , (9) 

and noise strength D cxt - The numerical integration is carried out by the stan- 
dard Euler algorithm with the step size At w 2 ■ 10~ 3 ms. The "Numerical 
Recipes" routine ran2 is used for the generation of independent random num- 



bers 21 1 with the Box-Muller algorithm providing the Gaussian distributed 
random numbers. The total integration time is chosen to be a multiple of the 
driving period Tq = 27r/f2, as to ensure that the spectral line of the driving 
signal is centered on a computed value of the power spectral densities. From 
the stochastic voltage signal V(t) we extract a point process of spike occurrences 

{**}: 

N 

u(t):=2*(t-t,), (10) 

8=1 

where N is the total number of spikes occurring during the elapsed time interval. 
The occurrence of a spike in the voltage signal V(t) is obtained by upward- 
crossing a certain detection threshold value Vq. It turns out that the threshold 
can be varied over a wide range with no effect on the resulting spike train 
dynamics. However, to account for the typical spike duration a time interval of 
2ms has been used. The power spectral density of the spike train (PSD U ), the 
interspike interval histogram (ISIH) and the coefficient of variation (CV) have 
been analyzed. The coefficient of variation CV, which presents a measure of the 
spike coherence, reads: 



where (T) := limjv^oo i ~ *») and (T 2 ) '■= lim^oo W XX^ ; +i ~ 

are the mean and mean-squared interspike intervals, respectively. From the 
PSD„ we obtain the height of the spectral line of the periodic stimulus as the 
difference between the peak value and its background offset. The signal-to-noise 
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Figure 1: The CV denned in equation ( [Ti"|) plotted versus area S (a) for 
A = 1.0/iA/cm 2 , = 0.3ms _1 , in the absence of external noise D ext = (solid 
line) and with no external stimulus applied (dotted line) . The ISIH are depicted 
in the presence of signal for S = 1/im 2 (b), and S = 32/im 2 (c). The vertical 
lines denote the driving period and the first two multiples. 



ratio (SNR) is then given by the ratio of signal peak height to the background 
height (in the units of spectral resolution of signals). 

We have analyzed the spike coherence in the autonomous, nondriven regime 
(i.e., we use Text = 0) as a function of the decreasing cluster size. Our simulation 
reveals, cf. fig. 0(a), the novel phenomenon of intrinsic coherence resonance , 
where the order in the spike sequence increases; i.e. the CV is decreasing, solely 
due to the presence of internal noise. The fully disordered sequence (which 
corresponds to a Poissonian spike train) would assume the value CV = 1. We 
note, however, that near S = 1/im 2 (optimal dose of internal noise), CV ~ 0.44, 
i.e. the spike train becomes distinctly more ordered! For S < 1/im 2 , the 
internal noise increases further beyond the optimal value and destroys the order 
in spiking again. It is worth noting that for S < 1/xm the model reaches limiting 
validity; in that regime the kinetic scheme Jn| [l^, [l9| should be used instead. 
Such a corresponding study, however, has been put forward independently by 
Jung and Shuai fll9| ; their results are in qualitative agreement with our findings. 
Next we switch on an external sinusoidal driving: Interestingly enough the 
interspike intervals distribution is not affected for small patch sizes, cf. fig. 0(b) 
for S — 1/im 2 . In this case, the spike-activity possesses an internal rhythm 
which dominates over the external disturbances. For larger patch sizes the 
internal noise decreases and the periodic drive induces a reduction of the CV 
as compared to the undriven case, note the right arrow in fig. 0(a). In this 
latter regime the external driving rules the spiking activity as depicted with the 
characteristic peaks in the ISIH in fig. 0(c) at multiple driving periods. 

Next, the focus is on the SNR in absence of external noise, see fig. 0(a). 
Here we discover the novel effect of genuine intrinsic stochastic resonance, where 
the response of the system to the external stimulus is optimized solely due to 
internal, omnipresent noise. For the given parameters it occurs at S w 32/im 2 . 
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Figure 2: The signal-to- noise ratio (SNR) for an external sinusoidal stimulus 
with amplitude A = 1.0/iA/cm 2 and frequency SI = 0.3ms _1 for different ob- 
servation areas: (a) No external noise is applied; (b) SNR versus the external 
noise for the system sizes indicated by the arrows in fig. ||(a). S = 8/zm 2 
(N K = 144, 7V Na = 480): solid line through the diamonds, S = 16/mr 2 
(JV"k = 288, A^Na = 960): long dashed line connecting the circles, S = 32/im 2 
(JVk = 576, A^Na = 1920): short dashed line through the triangles, S — 64/im 2 
(JVk = 1152, ATuja = 3840): dotted line connecting the squares. The situation 
with no internal noise (i.e., formally S — > oo) is depicted by the dotted line. 



For S < 32/im 2 growing internal noise monotonically deteriorates the system 
response. Under such circumstances, one would predict that the addition of 
an external noise (which corresponds to the conventional situation in biological 
SR studies) cannot improve SNR further, i.e. conventional SR will not be 
exhibited. Our numerical simulations, fig. ^|(b), fully confirms this prediction. 
Conventional stochastic resonance therefore occurs only for large membrane 
patches beyond optimal size, and reaches saturation in the limit S — > oo (limit 
of the deterministic Hodgkin-Huxley model). Thus, the observed biological SR 

H is rooted in the collective properties of large ion channels arrays, where 
ion channels are globally coupled via the common membrane potential V(t). 

In conclusion, we have investigated stochastic and coherence resonance in a 
noisy generalization of the Hodgkin-Huxley model for excitable biological cell 
membrane patches. The spontaneous fluctuations of the membrane conductiv- 
ity due to the individual channel dynamics has systematically been taken into 
account. We have shown that the excitable membrane patches with an obser- 
vation area around S ~ 1/im 2 exhibit a rhythmic spiking activity optimized by 
everpresent internal noise. In other words, the collective dynamics of globally 
coupled ion channels become more ordered solely due to internal noise. This 
new effect can be regarded as the intrinsic coherence resonance phenomenon; 
it presents a first important result of our work. This very finding has also been 
confirmed independently within a different approach by Jung and Shuai |jl9f 
(see the accompanying paper). A second main result of this study refers to the 
phenomenon of intrinsic SR; thereby the internal noise alone gives rise to a SR 
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behavior, see fig. ||(a). Conventional SR versus the external noise intensity takes 
place as well for sufficiently large membrane patches where the internal noise 
strength alone is not yet at its optimal value. We hence conclude that observed 
biological SR likely is rooted in the collective properties of globally coupled ion 
channel assemblies. 
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Peter Jung. 



References 

[I] Gammaitoni L., Hanggi P., Jung P. and Marchesoni F., Rev. Mod. Phys. 
70, 223 (1998). 

[2] Douglass J. K., Wilkens L., Pantazelou E. and Moss F., Nature (London) 
365, 337 (1993). 

[3] Levin J. E. and Miller J. P., Nature (London) 380, 165 (1996). 

[4] Collins J. J., Imhoff T. T. and Grigg P., Nature (London) 383 770 (1996). 

[5] Simonotto E., Riani M., Seife C, Roberts M., Twitty J. and Moss F., Phys. 
Rev. Lett. 78, 1186 (1997). 

[6] Hidaka L, Nozaki D. and Yamamoto Y., Phys. Rev. Lett. 85, 3740 (2000). 

[7] Bezrukov S. M. and Vodyanoy I., Nature (London) 378, 362 (1995); 385, 
319 (1997). 

[8] Goychuk I. and Hanggi P., Phys. Rev. E 61, 4272 (2000). 

[9] Petracchi D., Pellegrini M., Pellcgrino M., Barbi M. and Moss F., Biophys. 
J. 66, 1844 (1994). 

[10] Hodgkin A. L. and Huxley A. F., J. Physiol. (London) 117, 500 (1952). 

[II] Longtin A., J. Stat. Phys. 70, 309 (1993); Wiesenfeld K., Pierson D., Pan- 
tazelou E., Dames C. and Moss F., Phys. Rev. Lett. 72, 2125 (1994); Collins 
J. J., Chow C. C, Capela A. C. and Imhoff T. T., Phys. Rev. E 54, 5575 
(1996); Lee S.-G. and Kim S., Phys. Rev. E 60, 826 (1999). 

[12] Pikovsky A. S. and Kurths J., Phys. Rev. Lett. 78, 775 (1997); Lee S.- 
C, Neiman A. and Kim S., Phys. Rev. E 57, 3292 (1998); Lindner B. and 
Schimansky-Geier L., Phys. Rev. E 61, 6103 (2000). 

[13] Hille B., Ionic Channels of Excitable Membranes, 2nd ed. (Sinauer Asso- 
ciates, Sunderland, MA, 1992). 



8 



[14] Nossal R. J. and Lecar H., Molecular and Cell Biophysics, (Addison- Wesley, 
Redwood City, 1991). 

[15] Lowen S. B., Liebovitch L. S. and White J. A., Phys. Rev. E 59, 5970 

(1999) . 

[16] White J. A., Rubinstein J. T. and Kay A. R., Trends Neurosci. 23, 131 

(2000) . 

[17] Lecar H. and Nossal R., Biophys. J. 11, 1068 (1971). 

[18] Clay J. R. and DeFelicc L. J., Biophys. J. 42, 151 (1983); Strassberg A. F. 
and DeFelicc L. J., Neural Comput. 5, 843 (1993); DeFelice L. J. and Isaac 
A., J. Stat. Phys. 70, 339 (1993). 

[19] Jung P. and Shuai J., "Optimal sizes of ion channel clusters", preprint, 

(2001) . 

[20] Fox R. F. and Lu Y., Phys. Rev. E 49, 3421 (1994). 

[21] Press W. H., Tcukolsky S. A., Vetterling W. T. and Flanncry B. P., Nu- 
merical Recipes in C, 2nd ed. (Cambridge University Press, 1992). 



9 



